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0\ \ ABSTRACT 

Several recent papers have made claims about the detection of an asymmetric distribu- 
tion of large scale power in the cosmic microwave background anisotropy as measured 
j^ ' by the WMAP satellite. In this paper, we investigate how the estimates of particular 

fT^ \ cosmological parameters vary when inferred from power spectra computed separately 

CNJ ' on the northern and southern hemispheres in three frames of reference: the galactic, 

>0 \ the ecliptic and the specific frame of reference which maximises the asymmetry be- 

tween the power spectra from the corresponding hemispheres. Such a study is intended 
to quantify the consistency of the observed spectral variations in the context of the 
inflationary-scenario-inspired models with which the data are commonly compared. 
We focus our investigation on the three specific parameters to which the analysis is 
most sensitive - the spectral index n, the amplitude of fluctuations A and the optical 
O \ depth r - and flnd interesting variations in their estimates as determined from differ- 

■j-J ' ent hemispheres. When using a gaussian prior on the spectral index n centred at n = 1 

^ \ with a flat prior on the optical depth t, the preferred value for the optical depth (de- 

rived in the reference frame of maximum asymmetry) in the north is r = with upper 
fimit T < 0.08, whereas in the south we find r = 0.24too7 (68% confidence level). The 
^^ ■ latter result is inconsistent with r = at the 2a level. The estimated optical depth of 

H I T — 0.17 on the (nearly) full sky found by the WMAP collaboration and confirmed 

independently here, could thus in large part originate in structure associated with 
the southern hemisphere. Furthermore, putting a prior on r, we find values of the 
spectral index which are inconsistent between the opposite hemispheres. The exact 
values depend on the prior on r. Our conclusions remain unaltered even when, on the 
basis of putative residual foreground contamination, the multipole range ^ = 2 — 7 
is excluded from the analysis. While our results should not be considered more than 
suggestive, the significance of the parameter differences in the two hemispheres being 
typically of order 2ct, if they are confirmed with the higher sensitivity WMAP 2-year 
data then it may be necessary to question the assumption of cosmological isotropy 
and the conceptual framework within which studies of the CMB anisotropy are made. 
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1 INTRODUCTION 

Several recent papers have determined from an analysis of the data from the WMAP satellite lIBennett et al. 20031 that the 
fluctuations in the cosmic microwave background (CMB) are asymmetrically distributed on the sky (Park 2003 Eriksc n et al. 2004al 
IVielva et al. 20031 lEriksen et al. 200 4b Hans en et al. 2004al |Copi et al. 2003| IHansen et al. 2004bl iLarson fc Wandeh 20041 
ICruz et al. 20041 [Land fc Magueijo'2004j . In particular, it was shown (|Eriksen et al. 2004al IHansen et al. 2004b)l that the 
power spectra as computed in the northern and southern hemispheres (in both the galactic and ecliptic frames of reference) 
differ significantly on large angular scales. Specifically, for multipoles between 5 and 40, the southern hemisphere spectrum 
showed evidence of excess power whereas the northern spectrum was remarkably lacking therein, with respect to the canon- 
ical full sky estimate ^. The asymmetry in the spectrum over these multipoles was found to be significant at the 2 — 3(7 
level compared to Gaussian realizations of the best fit WMAP power spectrum. Moreover, the axis of maximum asymmetry 
was found to have its north pole at (80°, 57°) (galactic co-latitude and longitude). This direction was subsequently con- 
firmed by (IHansen et al. 2004al [Land fc Magueijo 2004| l using tests of non-Gaussianity. A similar (albeit at lower statistical 
significance) asymmetric distribution of power was also found using the data of the COBE-T)MK Ipcnnctt et al. 1996) exper- 
iment. There have also been detections of non-Gaussianity in the WMAP data jNaselsky et al. 2003, ,Naselsky et al. 2004a 



INaselsky et al. 2004b| IPark 20031 IVielva et al. 20031 iHaiisen et al. 2004al ILarson fc Wandelt"2004l [Land fc Magueijo 2004 



lEriksen et al. 2004bl ICruz et al. 2004^ . Whilst some of these have been claimed to be consistent with residual foreground 
nals ( |Naselsky et al. 2003| [Naselsky et al. 2004a| [Naselsky et al. 2004bt , it remains uncertain whether this non-Gaussianity 
is related to the power spectrum asymmetry. Indeed, although indications of a possible foreground contamination of the 
quadrupole and octupole have been found (i.e. lIHansen et al. 2004bllSchwarz et al. 2004llEriksen et al. 2004cl|Slosar fc Seljak 2004] !), 
the asymmetry at higher multipoles is harder to explain in terms of incorrectly subtracted foregrounds since it was found to 
be present at the same amplitude in all the WMAP frequencies and does not seem to decrease significantly with an increasing 
galactic cut IjEriksen et al. 2004al IHansen et al. 2004bfl . 

In this paper, the aim is not to obtain high-precision estimates of cosmological parameters in themselves, but rather to 
assess which parameters are more likely to be affected by the measured asymmetry and then to investigate and quantify the 
differences of the parameters between the north and the south. We will restrict the class of models for which we determine 
these parameterisations to the basic (inflationary-scenario-inspired) cosmological model as considered by the WMAP team 
l |Spergel et al. 2003| l - a flat ACDM model with power-law adiabatic primordial perturbations. The parameters of this model 
are the amplitude A and the spectral index n of the primordial power spectrum of scalar perturbations (modelled as Ak"), the 
optical depth r, the physical density of baryons and cold dark matter, Qbh^ and Q.cdmh? (where h is the Hubble constant in 
units of 100 km/s/Mpc), and the cosmological constant Q,a- Note that when fixing = 1, the Hubble constant h is set when 
the parameters Q,bh? , Qcdmh^ and Qa are given. We determine estimates of these parameters in the northern and southern 
hemisphere in three reference frames: galactic, ecliptic and that frame which maximises the asymmetry in the power spectra 
computed in the corresponding hemispheres. 



2 METHOD 

In lIHansen et al. 2004bL the power spectrum was estimated on differently oriented hemispheres using the Gabor transform 
formalism introduced by lIHansen et al. 20021 IHansen fc Gorski 2003t and extended in lIHansen et al. 2004bll . A brief outline 
of the method is provided here. 

A Gaussian likelihood ansatz is adopted with the pseudo power spectrum as input. The pseudo power spectrum Ce is 
given by 

m=-e 
where 

aim = / T{n)G{n)dh (2) 



are the spherical harmonic coefficients on the observed sky. Here, n is the position on the sky, T{n) is the CMB temperature 
fluctuation amplitude at that position and G{n) is the mask. The latter quantity takes zero value in those parts of the sky 
which are either unobserved or have been excluded from the analysis, and is otherwise one. The likelihood can now be written 
as 

^ This was considered to be the best fit power-law spectrum, as determined from the WMAP analysis using the Kp2 sky cut including 
the point source mask. Whenever we refer to a 'full-sky' analysis this coverage is implied, and the nomenclature is simply for convenience 
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„-idTM-^d 

C= ' (3) 

V27rdetM 

where the elements of the data vector are 

d, = Ct^-{Cl)-{C?^). (4) 

Here, the first term is the observed pseudo power spectrum including noise and the second term is the expected signal 
pseudo spectrum as a function of the underlying full sky power spectrum Ci. The last term contains the expected noise 
pseudo spectrum for the observed part of the sky. Similarly, the matrix expressing the correlation between the pseudo power 
spectrum coefficients can be written as 

M,j =< Cefii^ >-<Ci^>< Ce^ >= Mfj + Mf^ + Mf^, (5) 

where there is one contribution from the signal (S), one from the noise (N) and one from signal- noise mixing (X). The exact 
form of these matrices, their dependence on the full sky power spectrum C< and details of their computation are discussed 
in llflansen et al. 20021 Iflansen fc Gorski 20031 IJiansen et al. 2004b( l. In these papers, the power spectrum was estimated by 
maximising this likelihood with respect to Ci. As the likelihood is not Gaussian with respect to the power spectrum, there is 
a small bias introduced in the estimate of the d l |Bond, Jaffe fc Knox 2000l l at the smallest multipoles, for which a correction 
was found in IJiansen et al. 2004bt . In this paper, we maximise the likelihood with respect to the cosmological parameters. We 
have tested this method on simulated maps with known cosmological parameters and have demonstrated that the procedure 
gives an unbiased estimate of the input cosmological parameters. 



3 PREPARING THE DATA 

In this paper, we analyse the publicly available WMAP data^ after correction for Galactic foregrounds. We consider data 
from the V and W frequency bands only since these are less likely to be contaminated by residual foregrounds. The maps 
from the different channels are coadded using inverse noise- weighting IIHinshaw et al. 2003t : 

where c runs over the frequency channels V1,V2,W1,W2,W3,W4. All analysis is performed using the WMAP Kp2 sky cut, 
including the point source mask. 



4 THE RESULTS 

Since we are interested in estimating the relative difference in parameters between the hemispheres rather than an exact 
evaluation of their values to high precision, the multipoles used in the analysis are limited to the ranges which are important 
for characterising the difference between the northern and the southern spectra. The asymmetry manifests itself as both a 
difference in shape over the first 40 multipoles and as a difference in amplitude between the low-order multipoles and those 
at higher l. Therefore, multipoles in the range ^ = 2 — 60 are selected to account for the shape of the spectrum at large 
angular scales, together with a set of multipoles centred around the first two peaks {I — 200 — 240 and £ — 520 — 560) 
to constrain the shape of the spectrum at smaller scales. As noted in the introduction, we are generally interested in the 
parameters A, n, r, Q,bh^, i^A and ftcdmh^ , but have found that the asymmetry does not depend strongly on the exact value 
of the latter two values. This should not be surprising, since when Q is fixed, the parameters f^A and Q,cdmh^ influence the 
position of the acoustic peaks in the power spectrum, which does not differ significantly between the different hemispheres 
IjHansen et al. 2004b^ . Consequently, it was decided for simplicity to fix these parameters to the best fit values found by the 
WMAP team (Qa = 0.74 and ftcdmh^ = 0.11). These values are supported by independent astronomical observations that 
are consistent with the WMAP findings (for example, large scale structure matter distribution llVerde et al. 20021 . HST Key 
Project, lIFreedman et al. 200111 . type la supernovae lIRiess et al. 200l1lPerlmutter et al. igogt ) and that are not likely to be 
affected by the asymmetry. In addition, we have explicitly checked that the results in this paper are robust with respect to 
changes in these parameters. 

For the analysis of the remaining four parameters, it was found that the adoption of a flat prior resulted in a poorly 
behaved likelihood surface, with a large number of false local maxima. This is a consequence of strong degeneracies between the 
parameters. However, by setting stronger priors on two of the parameters, these degeneracies were reduced and we were able 
to obtain estimates of the remaining two parameters. For the baryon fraction we adopted a Gaussian prior with Oj/i^ — 0.022 

^ These can be obtained at the LAMBDA website: http://lambda.gsfc.nasa.gov/ 
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and a — 0.001. This choice was motivated by observations of the primordial deuterium abundances which give values in the 
range Q.^h'^ = 0.020 - 0.025 ID'Odorico et al. 200l1lO'Meara et al. 2001IIPettini fc Bowen 200111 . Ahernative Gaussian priors 
centred at different values in this range and with different a do not qualitatively affect the results presented below. Finally, 
two sets of parameter estimates were computed: 

• Gaussian prior on the spectral index n. Two different priors were assumed, one centred at n = 1 corresponding to 
a flat Harrison- Zeldovich spectrum, the second at n = 0.96 as suggested by the best fit WMAP spectral index, both with a 
dispersion a — 0.03. We found the results robust to changes in a. t and A were then estimated with a fiat prior on these two 
parameters. The results of this procedure as applied to the 6 different hemispheres are shown in table (Q and Q and figure 

• Gaussian prior on the optical depth r. Again two priors were adopted, one centred at r = and one centred at 
the WMAP best fit value r = 0.17, both with o = 0.05. The qualitative results and conclusions were found not to change 
significantly under changes in o. The results of the parameter estimation procedure on the 6 different hemispheres are shown 
in tables and @ and figure ^. 

The estimated parameters given in the tables are found in the following way: 

(i) For a given parameter, we found the marginalized likelihood by integrating over all other parameters. 

(ii) The best fit value quoted in the tables is the value which maximises the marginalized likelihood. This value was found 
to be consistent with the value which maximises the full likelihood also given in the tables. 

(iii) The la and 2a ranges given in the tables are those for which the integrated marginalized likelihood is 68% and 95% of 
the total integral respectively. Note that because of degeneracies between parameters, some of these error bars are asymmetric. 
Forcing a stronger prior (lower a) on the other parameters, the error bars become more symmetric. 

. The figures were created with following procedure: 

(i) For each possible pair of two parameters we determined the likelihood marginalized over the remaining parameters. 

(ii) The contours shown in the plot are the areas for which the integral of the two dimensional marginalized likelihood is 
less than 68% of the total. The shaded zone shows the full sky results, the zones filled with horizontal lines with a dashed 
border show the results on the northern hemisphere and the zones filled with vertical lines with dotted borders show the 
results on the southern hemisphere. 

(iii) The position of the maximum of this two-dimensional likelihood is shown as a filled square (full sky), open square 
(northern hemisphere) and cross (southern hemisphere). 

The most striking feature in these tables is the fact that when assuming a Harrison-Zeldovich spectrum with a spectral 
index close to one, the optical depth in the south is found to be r ~ 0.23 whereas the most likely optical depth in the north is 
consistent with zero. This holds for all three frames of reference. If instead we use a flat prior on n and put a Gaussian prior on 
r, the value of the spectral index becomes significantly different in the northern and the southern hemispheres with a larger 
value in the south, the exact value depending on where the prior on r is centred. Note further that, for the full sky, the estimate 
of the optical depth (r = 0.16 assuming a spectral index close to 1) determined from the marginalized likelihood is close to 
the value found by the WMAP team, despite the fact that the WMAP results were derived additionally using polarisation 
data which is more sensitive to the optical depth than the pure temperature data used here. It should be noted, however, 
that as a consequence of strong degeneracies, the marginalized likelihood in r is very flat and the maximum at r = 0.16 is 
not a very peaked maximum. This is also seen in figure^where the la contour is very extended in the r direction. This also 
explains why there is such a large difference between the estimated value of r = 0.16 using the marginalized likelihood and 
the estimate of r = using the full likelihood as seen in table Q . As will be discussed below, when excluding the very first 
multipoles the preferred value of r on the full sky increases. 

We will now outline how the results obtained above can be understood qualitatively by considering the manner in which 
the parameters of interest. A, r and n infiuence the shape of the spectrum. First, we will consider the case when n is fixed 
and discuss the variation of A and r. The optical depth r attenuates the amplitude of the higher multipoles (from I ~ 180 
and above) by a factor of roughly e~^'^. The amplitude and the optical depth cannot then be estimated separately over this 
range, so an effective resulting amplitude is written as Ae~'^'^ . We know that the power spectra on these scales are consistent 
between the different hemispheres, and therefore the estimates of A and r should lie close to the line Ae"^'" — constant in the 
A — T diagram, as is indeed seen in figure Q. In order to break this degeneracy, the shape and amplitude of the lower order 
multipoles are needed, since the optical depth has only a small infiuence on the amplitudes in this range. Therefore, when 
fixing the spectrum at higher multipoles, a greater amplitude for the low multipoles with respect to the higher ones yields a 
higher optical depth and vice versa. In IIHansen et al. 2004b|l . it was shown that the amplitude of the spectrum for the low 
multipoles {£ — 2 — 40) is higher in the south than in the north, one could thus expect a higher value of r in the south. This 
is exactly what we have seen in tables 11I2II and in figure (0. 

Now we consider the shape of the spectrum in the first multipoles, where a strong degeneracy between r and n manifests 
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Table 1. Estimated parameters with a prior on n centred at n = 1.0 

Hemisphere optieal depth t Ict range (t) 2ct range (x) abs. max (7-) amplitude A \ct range (A) 2ct range (A) abs. max (A) 



Galactic north 


0.01 


0.0-0.15 


0.0-0.24 


0.03 


0.66 


0.60-0.85 


0.57-1.03 


0.65 


Galactic south 


0.23 


0.16-0.28 


0.07-0.31 


0.21 


0.97 


0.85-1.09 


0.69-1.15 


0.95 


Ecliptic north 


0.0 


0.0-0.11 


0.0-0.21 


0.03 


0.66 


0.60-0.78 


0.57-0.95 


0.65 


Ecliptic south 


0.22 


0.15-0.28 


0.06-0.31 


0.18 


1.06 


0.92-1.12 


0.74-1.27 


1.00 


Max. ass. north 


0.0 


0.0-0. OS 


0.0-0. IS 


0.03 


0.65 


0.59-0.74 


0.56-0.88 


0.65 


Max. ass. south 


0.24 


0.17-0.3 


0.08-0.32 


0.24 


1.03 


0.89-1.12 


0.71-1.18 


1.05 


Full sky 


0.16 


0.01-0.22 


0.0-0.26 


0.0 


0.67 


0.64-0.93 


0.61-1.10 


0.65 



Table 2. Estimated parameters with a prior on n centred at n = 0.96 

Hemisphere optical depth t Ict range (x) 2ct range (x) abs. max (x) amplitude A Ict range (A) 2ct range (A) abs. max (A) 



Galactic north 


0.0 


0.0-0.1 


0.0-0.2 


0.0 


0.62 


0.58-0.73 


0.66-0.89 


0.60 


Galactic south 


0.17 


0.08-0.23 


0.0-0.26 


0.15 


0.82 


0.68-0.95 


0.58-1.04 


0.80 


Ecliptic north 


0.0 


0.0-0.07 


0.0-0.17 


0.0 


0.62 


0.58-0.71 


0.66-0.83 


0.60 


Ecliptic south 


0.15 


0.06-0.22 


0.0-0.25 


0.00 


0.88 


0.71-1.01 


0.64-1.16 


0.80 


Max. ass. north 


0.0 


0.0-0.06 


0.0-0.14 


0.0 


0.61 


0.58-0.68 


0.66-0.79 


0.60 


Max. ass. south 


O.IS 


0.09-0.24 


0.0-0.27 


0.18 


0.86 


0.70-0.98 


0.59-1.08 


0.90 


Full sky 


0.0 


0.0-0.12 


0.0-0.21 


0.03 


0.66 


0.60-0.79 


0.57-0.96 


0.65 



itself. Specifically, when the optical depth is increased, the shape of the spectrum is tilted in such a way that the lowest I 
part of the spectrum is raised with respect to the rest of the low I plateau. Lowering the spectral index n results in similar 
behaviour. As shown in ( Hansen et al. 2004b ), the southern spectra have this tilt favouring either lower n or higher r values. 
Both of these effects are seen in the results presented here - the estimates of r with a prior centred at n = 0.96 are lower 
than the corresponding values determined with the prior centred at n = 1. Moreover, the estimates of the spectral index n 
axe higher when the prior on r is centred at r = 0.17 than when centred at r = 0. 

In order to understand the significance of the results further, 512 simulations were generated with the effective beam and 
noise of the co-added W-l-V channel, and with the input power spectrum specified by the WMAP best fit full sky parameters 
for a power law index. From these maps, it was determined that 5% showed evidence for a difference in r between the northern 
and southern hemispheres (defined in the maximum asymmetry reference frame of the WMAP data) as large as that observed 
here. However, only one of these 512 maps had a difference in the power spectrum comparable to that in the WMAP data. We 
then considered the subset of simulated maps which had, relative to the WMAP results, a larger difference in r but smaller 
difference in power spectra between the northern and the southern hemispheres. In most cases, the large differences in r were 
due to a high value in one hemisphere driven by a strong upward fluctuation in a few multipoles. This is qualitatively different 
from the WMAP behaviour, in which the spectrum on one of the hemispheres is systematically higher over a large range of I. 

Finally, we note that IIHansen et al. 2004bt found some evidence for possible residual galactic foreground contamination 
in the multipole range £ = 2 — 4. Therefore, the parameters of interest in this work were re-evaluated after excluding the 
multipoles 1 — 2 — 7 (this range was appropriate because of the binning used in this analysis). There has been considerable 
debate in the literature as to whether this multipole range is of abnormally low amplitude relative to the rest of the spectrum 
(and more particularly relative to the claimed best fit cosmological model). We do not add to this discussion here, but recall 
that such behaviour favours a lower value of r or spectral index n, so that on exclusion of these low-^ multipoles the preferred 
value of r and n should increase. This is indeed what is observed. However, the significance of the asymmetry and the difference 
in parameters between the north and the south change only slightly and the results are not reported here. 



Table 3. Estimated parameters with a prior on r centred at t = 0.0 

Hemisphere spectral index n Ict range {n) 2ct range {n) abs. max (ti) amplitude A Ict range (A) 2ct range (A) abs. max (A) 



Galactic north 


0.98 


0.95-1.01 


0.92-1.04 


0.97 


0.62 


0.58-0.68 


0.66-0.76 


0.60 


Galactic south 


0.93 


0.90-0.96 


0.87-0.98 


0.93 


0.61 


0.57-0.67 


0.64-0.76 


0.60 


Ecliptic north 


0.99 


0.97-1.03 


0.94-1.06 


0.98 


0.64 


0.59-0.70 


0.66-0.77 


0.60 


Ecliptic south 


0.96 


0.95-1.01 


0.92-1.04 


0.96 


0.71 


0.65-0.79 


0.62-0.89 


0.70 


Max. ass. north 


1.01 


0.99-1.06 


0.96-1.08 


0.99 


0.66 


0.60-0.70 


0.66-0.77 


0.60 


Max. ass. south 


0.93 


0.91-0.96 


0.87-0.99 


0.93 


0.62 


0.58-0.69 


0.66-0.77 


0.60 


Full sky 


0.96 


0.95-1.01 


0.92-1.04 


0.96 


0.65 


0.60-0.71 


0.57-0.79 


0.60 
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Figure 1. The 68% confidence levels for the marginalized likelihood determined using a Gaussian prior on the spectral index n with 
a = 0.03 and centred at either n = 1.0 (left column) or n = 0.96 (right column). The shaded area (solid border line) shows the interval 
for the full sky estimate, the area filled with horizontal lines (dotted border) shows the interval for the northern hemisphere and the area 
filled with vertical lines shows the interval for the southern hemisphere (dashed border). The filled square, open square and cross show 
the maximum for the full sky, northern hemisphere and southern hemisphere respectively. The reference frames are: galactic (upper row), 
ecliptic (middle row) and maximum asymmetry (lower row). Note the increasing disparity between the values derived from the northern 
and southern hemispheres as this sequence is descended. 
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Figure 2. The 68% confidence levels for the marginalized likelihood determined using a Gaussian prior on the optical depth r with 
a = 0.05 centred either at t = 0.0 (left column) or t = 0.17 (right column). The shaded area (solid border line) shows the interval for 
the full sky estimate, the area filled with horizontal lines (dotted border line) shows the interval for the northern hemisphere and the 
area filled with vertical lines shows the interval for the southern hemisphere (dashed border line). The filled square, open square and 
cross show the maximum for the full sky, northern hemisphere and southern hemisphere respectively. The reference frames are: galactic 
(upper row), ecliptic (middle row) and maximum asymmetry (lower row). 
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Table 4. Estimated parameters with a prior on r centred at r = 0.17 

Hemisphere spcetral index n Icr range (n) 2ct range (ti) abs. max (n) amplitude A Icr range (A) 2ct range (A) 



Galactic north 


1.02 


0.99-1.07 


0.95-1.11 


1.03 


0.89 


0.78-1.02 


0.69-1.15 


0.90 


Galactic south 


0.99 


0.94-1.02 


0.91-1.06 


0.99 


0.88 


0.77-1.00 


0.69-1.13 


0.90 


Ecliptic north 


1.06 


1.01-1.09 


0.97-1.13 


1.03 


0.89 


0.78-1.02 


0.69-1.12 


0.85 


Ecliptic south 


1.05 


1.01-1.08 


0.97-1.12 


1.05 


1.10 


0.98-1.25 


0.86-1.36 


1.10 


Max. ass. north 


l.OS 


1.03-1.11 


0.99-1.15 


1.08 


0.90 


0.79-1.02 


0.70-1.16 


0.95 


Max. ass. south 


D.98 


0.95-1.03 


0.92-1.07 


0.98 


0.91 


0.80-1.04 


0.72-1.17 


0.90 


Full sky 


1.04 


1.00-1.08 


0.96-1.12 


1.04 


0.94 


0.83-1.08 


0.74-1.21 


0.95 



5 DISCUSSION AND CONCLUSIONS 

Previous work has determined the presence of a significant asymmetric distribution of large scale power in the WMAP 
data. Such a result could be suggestive of incompletely understood foreground residuals or systematic artifacts, interesting 
physical effects taking place at post-recombination, or a more profound lack of the correct cosmological framework within 
which to understand the data. In an attempt to gain further insight, we have investigated the influence of the asymmetry on 
the inferred values of specific cosmological parameters computed in the northern and southern hemispheres of three frames 
of reference - the galactic, the ecliptic and the frame of reference which maximises the observed asymmetry of the power 
spectrum. It was determined that the parameters which are sensitive to variations in the shape of the power spectrum seen 
in the different hemispheres are: the amplitude of fluctuations A, the spectral index n and the optical depth r. However, 
significant degeneracies between the parameters were present, and in detailed work we elected to place prior constraints one 
one parameter, then consider the behaviour of the remaining two. 

After placing a Gaussian prior on n centred at n = 1 with a = 0.03, we found that the optical depth in the northern 
hemisphere was small (r < 0.08 at la confidence in the maximally asymmetric frame of reference), whereas in the south the 
values were of order r ~ 0.23, inconsistent with r = at the 2a level. Therefore, it seems possible that the main contribution 
to the optical depth r = 0.17 estimated on the full sky by the WMAP team (and independently re-estimated here to be 
r = 0.16) comes from the southern hemisphere. The WMAP value for the optical depth also includes information from the 
temperature-polarisation, T x E, cross-power spectrum. A crucial test of our hypothesis is then to make a similar hemisphere- 
partitioned analysis of the polarised signal. A varying optical depth in different directions might occur if the reionization of the 
universe were highly inhomogeneous, but it is not clear how to produce such a scenario. An alternate possibility might arise 
from the presence of a huge gas mass in the southern hemisphere. This possibility was also discussed by l)Cruz et al. 2004^ in 
the context of a strong non-Gaussian cold spot located in this region. However, the asymmetry in the spectrum remained even 
when this spot were masked out. Moreover, to explain such a huge difference in r, an enormous amount of gas distributed 
over a huge area would be required: assuming that a cluster of galaxies may have an optical depth of order r ~ 0.01 implies 
that the gas required to produce an amplitude of r ~ 0.23 would correspond to the content of several clusters along the line of 
sight in all directions on the southern hemisphere. If on the other hand, one assumes that the WMAP data are contaminated 
by some unknown source, then the surprisingly high value of r — 0.23 might be taken as an indication that the southern 
hemisphere is suspect. 

Further, the estimate of the amplitude of fluctuations A and the optical depth r using only the temperature power 
spectrum, are highly correlated and closely follow the line Ae~'^^ = constant. Consequently, a low estimate of r implies a low 
estimate of A, and vice versa. In the north where the optical depth is low, the amplitude is also found to be low {A ~ 0.7) 
whereas in the south the amplitude is found to be signiflcantly higher A ~ 1.0. As outlined in l |Spergel et al. 2003| l, some 
weak lensing experiments and studies using galaxy velocity flelds flnd values of the amplitude consistent with our northern 
hemisphere findings. This is also the case for a recent estimate of the amplitude using the cluster mass function obtained 
from the baryon mass of clusters of galaxies ()Voevodkin fc Vikhlinin 2003^ . However, other experiments such as studies of 
high redshift clusters and measurements of the SZ effect through the high £ tail of the GMB power spectrum (again see 
l |Spergel et al. 2003| l as well as a recent paper by (lOdman et al. 200411 ) indicate a higher value for the amplitude, consistent 
with our estimate on the southern hemisphere. 

Alternatively, if we put a Gaussian prior on r and adopt a flat prior on the spectral index, we find that the spectral 
index determined in the northern hemisphere is systematically higher than in the southern hemisphere, where the exact value 
depends on where the prior for r is centred. For a low value of r, we estimate low values inconsistent with n = 1 in the 
southern hemisphere. On the other hand, if r = 0.17 we find high values inconsistent with n = 1 in the northern hemisphere. 
If we consider that this behaviour is an indication that one of the hemispheres contains some unknown contamination, one 
cannot use this result to draw a conclusion about its location. 

Finally, we note that should the asymmetry remain unexplained in terms of residual or unknown foreground contamination 
or systematic artifacts, the possibility remains that it may be necessary to invoke new post-recombination physical processes 
or to accept that our current conceptual framework based on the principle of cosmological isotropy and inflationary-scenario- 
inspired models is inadequate. We await the release of the WMAP 2-year results, and particularly the polarisation information, 
to further test these possibilities. 



Cosmological parameters and the WMAP data revisited 9 



ACKNOWLEDGEMENTS 



We are grateful to D. Marinucci, P. Mazzotta, A. Vikhlinin and S.D.M. White for useful discussions and comments. FKH 
acknowledges financial support from the CMBNET Research Training Network. We acknowledge use of the HEALPix 
IjGorski et al. 1998^ software and analysis package and the CMBfast software package ( |Seljak fc Zaldarriaga 19 96') for de- 
riving the results in this paper. This research used resources of the National Energy Research Scientific Computing Center, 
which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC03-76SF00098. We 
acknowledge the use of the Legacy Archive for Microwave Background Data Analysis (LAMBDA) . Support for LAMBDA is 
provided by the NASA Office of Space Science. 

REFERENCES 

Bennett C. L. et al., 1996, ApJ Letters, 464, LI 
Bennett C. L. et al., 2003, ApJS, 148, 1 

Bond J. R., Jaffe A. H. & Knox L., 2000, ApJ, 533, 19 

Copl C. J., Huterer D. & Starkman G. D., 2003, preprint jastro-ph/031051H 

Cruz M., Martinez-Gonzalez E., Vielva P. & Cayon L., 2004, preprint jastro-ph/040534H 
D'Odorico S. et al., 2001, A&A, 368, L21 

Eriksen H. K., Hansen F. K., Banday A. J., Gorski K. M. & Lilje P. B., 2004a, ApJ, 605, 14 

Eriksen H. K., Novikov D. L, Lilje P. B., Banday A. J. & Gorski K. M., 2004b, preprint I a stro-ph/0401276t 
Eriksen H. K., Banday A. J., Gorski K. M. & Lilje P. B., 2004c, preprint I astro-ph/0403098j 
Freedman W. L. et al., 2001, ApJ, 553, 47 

Gorski, K. M., Hivon, E., Wandelt, B. D., in 'Analysis Issues for Large CMB Data Sets', 1998, eds A. J. Banday, R. K. Sheth and L. Da 
Costa, ESQ, Printpartners Ipskamp, NL, pp. 37-42 Ustro-ph/9812350j ; Healpix HOMEPAGE: , http://www.eso.org/science/healpix/ J 
Hansen F. K., Gorski K. M. & Hivon E., 2002, MNRAS, 336, 1304 

Hansen F. K. & Gorski K. M., 2003, MNRAS, 343, 559 

Hansen F. K., Cabella P., Marinucci D. & Vittorio N., 2004a, ApJL, in press, preprin t jastro-ph/03105li) 
Hansen F. K., Banday A. J. & Gorski K. M., 2004b, submitted to MNRAS, prcprint| astro-ph/0404206t 

Hinshaw G. et al., 2003, ApJSS, 148, 135 

Land K. & Magueijo J., 2004, preprint lastro-ph/0405519 I 

Larson D. L. & Wandelt B. D., 2004, preprint ( astro-ph/0404037 1 

Naselsky P. D., Doroslikcvich A. G. & Vorkhodanov O. V., 2003, ApJ, 599, L53 

Naselsky P. D., Doroshkevich A. G. & Vorkhodanov O. V., 2004a, MNRAS, 349, 695 

Naselsky P. D., Verkhodanov O. V., Chiang L. Y. & Novikov I. D., 20 04b, preprint jastro-ph/040552'3) 

Odman C. J., Hobson M., Lasenby A. & Melchiorri A., 2004, preprint jastro-ph/0405118t 

O'Mcara J. M. et al., 2001, ApJ, 552, 718 

Park C. G., 2004, MNRAS, 349, 313 

Perlmutter S. et al., 1999, ApJ, 517, 565 

Pettini M. & Bowen D. V., 2001, ApJ, 560, 41 

Riess A. G. et al., 2001, ApJ, 560, 49 

Seljak U. & Zaldarriaga M., 1996, ApJ, 469,437 

Schwarz D. J., Starkman G. D., Hute rer D. fc Copi C. J ., 2004, preprint jastro-ph/0403353t 
Slosar A. & Seljak U., 2004, preprint jastro-ph/0404567t 
Smoot G. F. et. al., 1992, ApJ Letters, 396, LI 
Spergel D. J. et al., 2003, ApJSS, 148, 175 

Verde L. et al., 2002, MNRAS, 335, 432 

Vielva P., Martinez-Gonzalez E., Barreiro R. B., Sanz J. L. & Cayon L., 2003, preprint ja5tro-ph/0310273^ 
Voevodkin A. fc Vikhlinin A., 2003, ApJ, 601, 610 



